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Abstract: The quantum mechanical mass of’t Hooft-Polyakov monopoles in the four¬ 
dimensional Georgi-Glashow is calculated non-perturbatively using lattice Monte Carlo 
simulations. This is done by imposing twisted boundary conditions that ensure there is one 
unit of magnetic charge on the lattice, and measuring the free energy difference between 
this ensemble and the vacuum. In the weak-coupling limit, the results can be used to 
determine the quantum correction to the classical mass, once renormalisation of couplings 
is taken properly into account. The methods can also be used to study the masses at strong 
coupling, i.e., near the critical point, where there are hints of a possible electric-magnetic 
duality. 
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1. Introduction 


’t Hooft-Polyakov monopoles ii are topological solitons in the Georgi-Glashow model Q 
and a wide range of other gauge field theories, including super Yang Mills theories and grand 
unified theories. They are non-linear objects in which energy is localised around a point 
in space and which therefore appear as point particles, and they carry non-zero magnetic 
charge. It is possible that these monopoles actually exist in nature, but so far they have 
not been discovered^ despite extensive searches |^. However, ’t Hooft-Polyakov monopoles 
are very important theoretically, because they provide a new way of looking at non-Abelian 
gauge field theories, complementary to the usual perturbative picture. In particular, this 
has shed more light on the puzzle of confinement [0, p. So far, concrete results have been 
limited to supersymmetric theories. 

The main reason for the lack of progress in non-supersymmetric theories is the diffi¬ 
culty of treating the quantum corrections to the classical monopole solution. For instance, 
calculating the quantum correction to a soliton mass is a complicated task. Even in simple 
one-dimensional models, it can typically only be calculated to one-loop order [P, and for 
’t Hooft-Polyakov monopoles the situation is even worse as only the leading logarithm is 
known |^]. This difficulty is avoided in supersymmetric models, because the symmetry 
protects the mass from quantum corrections. 

In this paper, the quantum mechanical mass of a’t Hooft-Polyakov monopole is cal¬ 
culated using lattice Monte Garlo simulations. The method was developed in Ref. |]R| 


^One candidate event Q was seen on Valentine’s Day 1982 but is unlikely to have been a real monopole. 
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and has been used earlier |]d] in a 2+1-dimensional model in which the monopoles are 
instanton-like space-time events rather than particle excitations. The mass is dehned using 
the free-energy difference between sectors with magnetic charges one and zero, and the 
corresponding ensembles are constructed using suitably twisted boundary conditions. This 
method has several advantages over the alternative approaches based on creation and anni¬ 
hilation operators [^, |^, ]ff| or fixed boundary conditions |^, ^|. In particular, it gives 


a unique, unambiguous result, since it requires neither gauge fixing, choice of a classical 
field configuration nor identification of individual monopoles in the field configurations. 

Analogous twisted boundary conditions have been used before to compute soliton 
masses in simpler models, such as 1+1-dimensional scalar field theory 0, 3+1-dimensional 
compact U(l) gauge theory |l^] and 2+1-dimensional Abelian Higgs model []T|. In the lat¬ 
ter case, the results provided evidence for an asymptotic duality near the critical point |p0|| : 
The model becomes equivalent to a scalar field theory with a global 0(2) symmetry, with 
vortices and scalar fields changing places. It is interesting to speculate whether an electric- 
magnetic duality might appear in the same way in the Georgi-Glashow model. These 
methods can, in principle, used to test that conjecture. 

The outline of this paper is the following: The Georgi-Glashow model and the classical 
’t Hooft-Polyakov solution are introduced in Section In Section the model is discretised 
on the lattice and the lattice magnetic field is defined. The twisted boundary conditions are 
discussed in Section In Sections ^ and ^ the classical and quantum mechanical monopole 
masses are computed, and the results are discussed in Section Finally, conclusions are 
presented in Section 0 


2. Georgi-Glashow model 

The 3+1-dimensional Georgi-Glashow model consists of an SU(2) gauge field and 
an Higgs field in the adjoint representation, with the Lagrangian 


1 , 


C = --TrF^^F^^^ + Tr[T)^, $][Z1^, d>] - m^Trd-^ - A (Trd*^) , 


( 2 . 1 ) 


where the covariant derivative and the field strength tension are defined as = 
+ igA^ and Dy\/ig. and ‘h are traceless, self-adjoint 2x2 matrices, they 

can be represented as linear combinations of Pauli a matrices. 


(Tl = 


0 1 
1 0 


0-2 = 


0 -i 
i 0 


0-3 = 


1 0 
0 -1 


( 2 . 2 ) 


as A^ = A“cr 


(J) = 

On classical level, the model has two dimensionless parameters, the coupling constants 
g and A, and the scale is set by m?. When m? is negative, the SU(2) symmetry is broken 
spontaneously to U(l) by a non-zero vacuum expectation value of the Higgs field Tr<I>^ = 
u^/2 = |m^|/2A. In the broken phase, the particle spectrum consists of a massless photon, 
electrically charged bosons with mass mw = gv, a neutral Higgs scalar with mass 
ruH = and massive magnetic monopoles |0, 0. The terms ’’electric” and ’’magnetic” 
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refer to the effective U(l) field strength tensor defined as 


- ^Tr$[Z)^, ^[D,, !>]. (2.3) 

In any smooth field configuration, the corresponding magnetic held Bi = eijkJ-jkl‘2' is 
sourceless (i.e., V • .6 = 0) whenever <h 7 ^ 0. This is easy to see in the unitary gauge, in 
which oc (T 3 , because Eq. ( |2.3| ) reduces to and therefore B = 'V x ^3. 

At zeros of <I>, the divergence is ± 47 r/gf times a delta function, indicating a magnetic charge 
of qm = dvr/^. 

The classical’t Hooft-Polyakov monopole solution |]i|, ^ is of the form 

= ^Higvr), 
gr^ 

Ai = -eaij^ [1 - K{gvr)] , (2.4) 


where H{x) and K{x) are functions that have to be determined numerically. It is easy 
to check that this solution is a magnetic charge in the above sense. Because the energy 
is localised around the origin, the solution describes a particle. Once the functions H{x) 
and K{x) have been found, it is easy to integrate the energy functional to calculate the 
mass of the particle, as it is simply given by the total energy of the configuration. The 
energy density falls as p ~ l/2g‘^r‘^, implying that the mass is finite but also that there is 
a long-range magnetic Coulomb force between monopoles, as expected. 

The classical monopole mass Md can be written as 


Mci = 


inmw 

5^ 


/(^) 


(2.5) 


where f{z) is a function oi z = mn/Tnw and is known to satisfy /(O) = 1 |^, ^|. It 


has recently been calculated numerically to a high accuracy |^|. Asymptotic expressions 
for small and large 2 ; had been found earlier \‘M, ^|, but the authors of Ref. reported 
that they had found an error in the small-z expansion. According to them, the correct 
expansion is 


f{z) = l + ]^z + ]^z 




( 2 . 6 ) 


For large z, they found that 


f{z) = 1.7866584240(2) - 2.228956(7)^"^ -h 7.U{l)z-^ + 0(z“3). (2.7) 


In quantum theory, the mass of a soliton can be defined as the difference between the 
ground state energies in the sectors with one and zero charge. In principle, it is possible 
to calculate this perturbatively to leading order Q. First, one needs to hnd the classical 
solution 4>o{x), and consider small fluctuations 5{t,x) around it, 

(l){t,x) = 4>o{x) + 6{t,x). ( 2 . 8 ) 
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When one drops higher-order terms in the Lagrangian, one is left with field 5 in an harmonic 
x-dependent potential U{6) = ^W'((/>o(x))5^. One needs to find the energy levels ujk of this 
field by solving the eigenvalue equation 


-v2 + y"(0o(x)) 


6k{x) = Loldkix). 


(2.9) 


The one-loop correction to the soliton mass is then given simply by the difference in the 
zero-point energies of one- and zero-soliton sectors, 

-^i-ioop = Mci + - ^ (o;^ — LOk) ) (2-10) 

k 


where coj. refers to the energies in the soliton background and in the trivial vacuum. 

One has to be careful with degeneracies and ultraviolet divergences, but the calculation 
can be carried out exactly in, for instance, the 1-|-1-dimensional model. In the presence 
of a kink, the energy spectrum consists of two discrete levels Wq = 0 and wf = (3/2)m^, 
and a continuum uj^ = {(f‘/2 -|- 2)m?. It is essential that one takes into account the same 
number of eigenvalues in the two sectors, and the best way to ensure is to do the calculation 
on a finite lattice and take the lattice spacing to zero and the lattice volume to infinity 
afterwards. This gives the result 


A^kink 


2V2 I f Vs 
^ T ^ yV 




( 2 . 11 ) 


The one-loop calculation of the monopole mass would go along the same lines, but 
there are many extra complications, which make it technically more difficult. Instead of 
one held, one has to consider two coupled helds. The background solution is not known 
analytically except in the special case of A = 0, and even then, the eigenvalue equation 
cannot be solved analytically. It is difficult to regularise the theory without breaking either 
gauge or rotation invariance, and in any case, the theory has much stronger ultraviolet 
divergences. Nevertheless, because the monopole mass is a physical quantity, it is hnite once 
the couplings have been renormalised. No separate renormalisation of the monopole mass 
is needed, and that means that the scale dependence of the resulting one-loop expression 
for the monopole mass would automatically be such that it cancels the running of the 
couplings. 

So far, only the leading logarithmic quantum correction near the BPS limit has been 
calculated ISI, 


M = 


47rmw 

5^ 


Stt^ 


m 


1 + Vj? ( — T' + ) + 0{gV 


m 


w 


( 2 . 12 ) 


An interesting aspect of this result is that the it is logarithmically divergent in the BPS 
limit. This is related to the Coleman-Weinberg effect | |26|| , which makes it impossible to 
actually reach the BPS limit in the quantum theory. Quantum corrections give rise to a 
logarithmic term V log 4’^ which means that if one tries to decrease the scalar mass below 
a certain value, the vacuum becomes unstable. This leads to a constraint mn > gmw- 
If one wants to be able to test Eq. (|2.12) numerically, the logarithmic term has to be 
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much larger than the constant term next to it, but that is only possible if g is very small. 
This, however, means that the whole quantum correction will be small and therefore more 
difficult to measure. Furthermore, having a large mass hierarchy such as mn ^ 'nnw 
means that one would need to use a very large lattice. For these reasons, a numerical test 
of Eq. (|2.12|) is not attempted in this paper. 


3. Lattice discretisation 


To study the model numerically, let us carry out a Wick rotation to Euclidean space and 
discretise the model in the standard way, 

Ce = 2'Y^ [Tr$(x)2 - Tt^{x)U^,{x)^{x + jl)Ul{x) 

+4 [2 - Tr[/^^(f)] + m^TV + a(TV (3.1) 

The scalar field $ is defined on lattice sites and the gauge field is represented by SU(2)- 
valued link variables [/^, which correspond roughly to The plaquette is 

defined as U^y{x) = U^{x)Uy{x + jl)ul{x + i>)ut{x). 

It will be crucial that the magnetic field can be defined on the lattice and that magnetic 
monopoles are therefore absolutely stable objects |^. This is highly non-trivial, because 
many other topological objects such as Yang-Mills instantons are not well defined on the 
lattice 1^ . To define the discretised version of the field strength tensor , note that the 
set of configurations with < 1 > = 0 at any lattice site is of measure zero and therefore these 
configurations do not contribute to any physical observables. One can therefore define a 
unit vector valued field ^ This expression makes sense because is always 

proportional to the 2 x 2 identity matrix 1. 

Because = 1, one can define a projection operator 11+ = (1 -|- ^)/2. Let us use it 
to define the projected link variable 


Uf,{x) = n+(x)[/^(x)n+(x -h ft), (3.2) 

which is essentially the compact Abelian gauge field that corresponds to the unbroken U(l) 
subgroup. The corresponding Abelian field strength tensor is 


afj,y = - aigTvu^{x)uy{x + gWJx + i>)ul(x), 


9 

and the lattice version of the magnetic field 




(3.3) 


(3.4) 


The lattice magnetic field Bi is a well-defined, gauge-invariant quantity. The magnetic 
charge density is given by its divergence, 

^ . r . . 1 ziTT 

(3.5) 


dvr 

G —Z, 
9 


Pm{x) = ^ \Bi{x + i)- Bi{x) 
i=l 

and note that it is quantised. Being a divergence of a vector field, the magnetic charge is 
automatically a conserved quantity. 
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4. Twisted boundary conditions 

Because the magnetic charge Qm = f d^xpM defined in the previous section is a well- 
defined, gauge-invariant, quantised and conserved quantity even on the lattice, it is a well- 
defined question to ask what the lowest energy eigenvalue with Qm = j9 is. Furthermore, 
since the total magnetic charge inside a volume is given by a surface integral over the 
boundary, one can fix the total charge in a simulation by choosing appropriate boundary 
conditions. 

In practice, one can therefore define separate partition functions for each magnetic 
charge sector, 

Zn= j F>C/;.Z)$exp(-5[17^,4>]), (4.1) 

Jn 

where the boundary conditions for each sector are such that they fix the magnetic charge to 
Qm = he., the net number of monopoles is N. Since monopoles of the same charge 

are not expected to form bound states, and since their interaction potential decreases with 
distance as 1/r, they will be non-interacting provided that the lattice is large enough. 
Denoting the length of the time direction by T, the partition function is therefore 

Z^ = exp(-|iV|MT)Zo, (4.2) 

where M is the quantum mechanical mass of the monopole, and Zq is the partition function 
with = 0. In particular, one can express the mass as 

M = -iln|l. (4.3) 

1 Zq 

It was shown in Ref. that this can be achieved by using suitably ’’twisted” boundary 
conditions. It is clear that periodic boundary conditions will not be useful, because they 
will fix the total charge to zero. On the other hand, they have the attractive feature 
that they preserve translation invariance and therefore, as all lattice points are equivalent, 
there will be no physical boundary. However, this does not require perfect periodicity: 
Periodicity up to the symmetries of the Lagrangian is enough. An obvious alternative is 
to use C-periodic boundary conditions |2^, 

U^{x -b Nf) = U*{x) = a 2 Uf,{x)a 2 , 

^{x + Nj) = <I>*(x) = —a2^(x)a2. (4.4) 

They flip the sign of the magnetic field as one goes through the boundary and therefore 
allow non-zero magnetic charges. 

In fact, it turns out that C-periodic boundary conditions allow the magnetic charge 
to have any even value [^. This means that the one-monopole sector is still excluded, 
and also that in practice the magnetic charge will always be zero, because as long as M is 
non-zero and T is chosen to be large enough (i.e., T ^ 1/M), 

OO 

^c= Z2fe = Zo (1 + 0(e-2^^)) . (4.5) 

k=—oo 
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When the boundary conditions in Eq. {^A) are written in the matrix form, it becomes 
obvious that one could have used cxi or <73 instead of 02 ■ They are all related to each other 
by gauge transformations and therefore describe identical physical situations. However, 
this observation allows one to define twisted boundary conditions 


Uf,{x + Nj) = crjUfj,{x)crj, 
^{x + Nj) = —aj^{x)aj, 


(4.6) 


which are locally equivalent to Eq. (4.4), but not globally. It is possible to carry out a 


gauge transform to turn the boundary conditions to Eq. (4.4) in any single direction, but it 
is not possible to do it to all three directions simultaneously. Considering the total charge 
of the lattice, one finds that these twisted boundary conditions only allow odd values 
and therefore 


Ztw = 


Z2k+i — Zi (^2 -h 0(( 


-2MT 


))• 


(4.7) 


k=—oo 


Thus, the ratio of the twisted and C-periodic boundary conditions can be used to calculate 
the monopole mass. 


1 In — = M - 
T Zc 


In 2 


+ 0(e-2^^) ^ M as T ^ 00 . 


(4.8) 


As such, this expression is of little use, because it is not possible measure partition 
functions directly in Monte Carlo simulations. One cannot write the ratio of partition 


functions in Eq. (4.8) as an expectation value either, because Ztw and Zq have different 
boundary conditions. One possible way to avoid this problem is to change the integration 
variables in Ztw in such a way that the new variables satisfy C-periodic boundary condi¬ 
tions. This changes the integrand, or equivalently the action S ^ S + AS. This way, one 
can express Eq. ( |4.8D in terms of an expectation value ZtwjZ^ = (exp(—AS))c, where the 
subscript C indicates that the expectation value is calculated with C-periodic boundary 
conditions. In principle, this is measurable in the simulations. The shift AS consists of 


line integrals of the magnetic field around the lattice |K]. In practice, this approach does 
not work, because exp(—AA) has very little overlap with the vacuum and one would need 
extremely high statistics to obtain any meaningful results. 

Let us, however, adopt a different strategy. Going back to Eq. ( |4.8D , we can differentiate 
the mass with respect to some parameter x. 


dM 

dx 


1 

T 


Aw \ / c 


(4.9) 


where the subscripts ”tw” and ”C” refer to expectation values calculated with twisted and 
C-periodic boundary conditions, respectively. If one start at a point where one knows the 
monopole mass, one can integrate this to obtain the mass at any other parameter values. 
Possible choices for the start point of the integration are the classical limit, where M can be 
calculated directly, or any point in the symmetric phase where the monopole mass vanishes. 
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Let us choose the latter option and use x = m?. Thus we can write 

dM 


dm? 




(4.10) 


If one chooses a large enough initial value for it is guaranteed to be in the symmetric 
phase. In fact, since one can only carry out a finite number of measurements, it is better 
to use finite differences instead of derivatives 


M(m 2 ) — M{m\) = — — In 




,tw 


-(m2-m2)X]^Tr#2\ 


(4.11) 


ml,C 


where the subscript m\ indicates that the expectation value is calculated at m? = m\. The 
spacing between different values of m? has to be fine enough so that the expectation values 
can be calculated reliably. 


5. Classical Mass 

It will be interesting to compare the measured quantum masses with classical results to 
determine the quantum correction. However, the quantum mass will be computed on a 
finite lattice, and therefore it does not make sense to compare it with the infinite-volume 
continuum expression (2.5). Instead, one needs to know the classical mass on the same finite 
lattice. That is straightforward to compute by minimising the lattice action 5tw ( |3.1|) with 
twisted boundary conditions. The C-periodic boundary conditions are compatible with 
the classical vacuum solution, and therefore the minimum action in that sector is simply 
5“™ = — (m^/4A)TL^. The classical mass is therefore given by 


Md = 


cmin cmin 
^tw 


omin ^4 

4A^ ■ 


(5.1) 


In fact, since the classical monopole configuration is time-independent, it is enough to have 
only one time step in the time direction, T = 1. 

The classical monopole mass was measured on three different lattices, 16^, 24^ and 
32^ using couplings A = 0.1 and g = l/\/5, which correspond to 2 ; = 1. The results are 
shown in Fig. [I|. The coloured solid lines correspond to different lattice sizes, the smallest 
being at the bottom. The top dashed line (black) is the infinite-volume mass given by 
Eq. ( |2.5|) with /(I) « 1.238 as computed in Ref. ||^. These results show a significant 
finite-size effect and demonstrate why it is necessary to compare the quantum result with 
the classical mass on the same lattice. 

Deep in the broken phase, where m? <C 0, the finite-size effect should be due to the 
magnetic Coulomb interaction between monopoles. Because our boundary conditions have 
the physical effect of charge conjugation, we effectively have monopoles and antimonopoles 
alternating in a cubic array, with distance L between them. The energy of such a configu¬ 
ration is 


E{L) = M - 


27r 




M - 


10.98 


(5.2) 









Figure 1: Classical monopole mass. 


The lower dashed lines (coloured) show the predicted finite-size effects for the relevant 
lattice volumes, and one can see that the agreement is good deep in the broken phase. In 
fact, the lattice values are slightly below the continuum results based on Ref. |23|. This is 
most likely due to discretisation effects. 

Although the Coulomb interaction describes the finite-size effects very well deep in 
the broken phase, it fails badly as m? —> 0. What happens there is that the size of the 
monopole, which is proportional to 1 /grows and eventually becomes comparable 
with the size of the lattice. At some point it becomes energetically favourable for the whole 
system to remain in the symmetric phase. Because the field is zero, the twisted action 
in Eq. vanishes, and the result is 


Esyrani{L) — V {0)L^ — —. (5-3) 

This is shown as a dotted line for L = 16 in Fig. |^, and agrees well with the result near 
= 0. At intermediate values of m^, the minimum energy configuration corresponds to a 
deformed monopole, and therefore the actual result interpolates smoothly between the two 
behaviours. Nevertheless, we have identified the main sources of finite-size effects in the 
classical calculation, and we are therefore in a position to compare quantum and classical 
calculations. 


6. Simulations 

In the quantum simulations, the ensembles of configurations with twisted and C-periodic 
boundary conditions were generated using the Metropolis algorithm. Three different lattice 
sizes were used: 16^, 24^ x 16 and 32^ x 16. The system was first equilibrated by carrying 
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out 20000-60000 update sweeps depending on the lattice size and the value of m?, and after 
that, measurements were carried out every 100 updates. The number of measurements for 
each value of vn? was between 100 and 1700. 



2 

m 


Figure 2: Quantum monopole mass (points) compared with the classical mass (lines). 


The expectation values needed for the change of M from m\ to can be calculated 
in two different ways, 




= ( e 




( 6 . 1 ) 


This was be used to check that the system was properly in equilibrium, the statistics were 
sufficient and the spacing between different values of was small enough. Dehning 


the change in the monopole mass (4.11) can be written as 


M{ml) - M{ml) = - (/gtw + / 2 ,tw - /i,c - / 2 ,c) • 


(6.3) 


The statistical error A/ in each fi^x was estimated using the bootstrap method and it was 
made sure that fi^x and f 2 ,x agreed within the errors. The error in the mass difference 
was estimated to be 


A [M(m2) — M(m^)]^ 


+^/i,c + ^/l.c + (/i,c — / 2 ,c)^ 


(6.4) 
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The differences were then summed up, starting from tuq, the highest value of m?, where 
M was assumed to be zero, 


N-l 

M{m%) = Y, - M{ml)) . (6.5) 

n=0 

The total error was calculated by assuming that the errors in the individual mass differences 
were independent. 


N-l 

A[M(m^)]2 = Y^ [M{mY) - M{ml)f (6.6) 

The results are shown in Fig. Note that the errors are highly correlated. In Fig. we 
show the derivative of the mass calculated from Eq. ( [1. 10 ). 


1500 r 


1000 


32 x16 
24^16 
16^ 


I 


500 


i 


» I • * f 


= s*f * 


* 




i 33; 
3: 


L 


-0.4 


-0.35 -0.3 

2 

m 


-0.25 


Figure 3: Derivative of the monopole mass. 


7. Discussion 

The mass derivative in Fig. |3| has a sharp peak, above which it drops rapidly to zero. This 
is compatible with the classical result dM / dm? oc 1/ \/—for negative m? and zero for 
positive values. As the horizontal lines show, the peak height is proportional to the linear 
size of the lattice, which is what one would expect to happen in a second-order phase 
transition. A fit to the peak position gives an infinite-volume value ~ 0.268 for the 
critical point. 

The curves in Fig. ^ show the classical results shifted by this amount. The quantum 
measurements agree fairly well with them near the critical point, but start to deviate 
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deeper in the broken phase. The qualitative behaviour, as well as the finite-size effects, are 
nevertheless similar. 

To really carry out a quantitative comparison of the quantum and classical results 
and to extract the quantum correction to the mass, one has to consider the renormali¬ 
sation of the parameters. The values of m?, A and g that were used in the simulations 
correspond to bare couplings, but one should compare the measurements with the classi¬ 
cal mass calculated using the corresponding renormalised couplings Xr and gn. The 
values of the renormalised couplings depend on the renormalisation scheme and scale, and 
therefore there is no unique way to compare the results. Furthermore, if one calculates 
the renormalisation counterterms to a certain order in perturbation theory, the value of 
quantum correction obtained by subtracting the classical value from the quantum result 
is only valid to the same order, even though the quantum mass itself has been calculated 
fully non-perturbatively. 



Figure 4: Measurements of the monopole mass (black circles) at = —0.35 with different lattice 
sizes L. The solid line is a fit of the form (5.2) and gives gn — 0.40(6). The blue crosses show the 
corresponding classical masses. 


It would, therefore, be best to use a physically meaningful renormalisation scheme and 
compute the renormalised couplings non-perturbatively. This can be done by choosing 
three observable quantities X, Y and Z, one for each coupling, and measuring their values 
(X), (Y) and {Z) in Monte Carlo simulations. One would then calculate the same quantities 
in the classical theory, and fix the values of the renormalised couplings by requiring that 
the classical values agree with the quantum measurements, 

Xci{m%,XR,gR) = {X) etc. (7.1) 

It would be natural to choose the masses of the perturbative excitations mn and mw 
as two of these observables, although measuring mw is non-trivial because of its electric 
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charge. One can choose the monopole charge as the third observable, because its value can 
be determined relatively straightforwardly from the finite-size effects of the monopole mass. 
In Fig. ^ we show the measured finite-size effects at m? = —0.35. A fit to Eq. ( ^.21) with g 
and M as free parameters gives qr = 0.40(6). It agrees with the bare value g = l/\/5 ~ 
0.447 within errors, so one would need better statistics to be actually able to measure the 
renormalisation counterterm. 

The change of the monopole mass as a function of is also directly related to the 
renormalisation of the theory. In a classical continuum theory, m? only fixes the scale, and 
dimensionless observables do not depend on it. In the quantum theory, this scale invariance 
is broken, and taking towards the critical point corresponds to renormalisation group 
flow towards infrared. Roughly speaking one can identify niH with the renormalisation scale 
/i. In principle, one should therefore be able to use the non-perturbative renormalisation 
scheme discussed above to follow the running of the couplings even in the non-perturbative 
regime near the critical point. 

One can speculate on what may happen based on the perturbative running of the 
couplings. The one-loop renormalisation group equations are 

~ ^‘^9r^R + Sgfl dg\ _ ^Or 

d log g, Svr^ ’ d log g Svr^ 

Moving towards infrared, \r decreases and g"^ increases. In fact, \r becomes negative at 
a non-zero g, i.e., before one reaches the critical point. This is a sign of the Coleman- 
Weinberg effect and means that there is a first-order phase transition. However, if gR has 
become large enough before this happens, the one-loop approximation is not valid any 
more, and it is possible that the critical point can be reached. This would mean that the 
line of first-order Coleman-Weinberg phase transitions ends at a tricritical point. Beyond 
that there is a second-order phase transition, around which the theory is strongly coupled. 

This is exactly what happens in the Abelian Higgs model in 2-|-l dimensions. There are 
strong arguments that in that case, the second-order phase transition has a dual description 
in terms of a global 0(2) model [^]. In the duality map, vortices and the fundamental 
scalar fields of the models change places. This was recently tested by measuring the critical 
behaviour of the vortex mass using a technique that was very similar to what has been 
discussed in this paper [p^]. The critical exponents that characterise that behaviour of 
vortices near the transition point were found to agree with the known critical exponents of 
the 0(2) model, which provides strong numerical evidence for the duality. 

If the Georgi-Glashow model has a second-order phase transition, it may have an 
analogous dual description. The one-loop renormalisation group equations suggest that as 
the critical point is approached, gR diverges and Xr goes to zero. The masses of the H 
and bosons and the monopoles should behave as 

1/2 1 

mn (X Xj{ , mw x gR, M (x g]^ , (7.3) 

implying that near the critical point, the bosons become much heavier than the other 
degrees of freedom and decouple. The Higgs scalar is neutral, and therefore it decouples 
as well. Thus one is left with massive magnetic monopoles coupled to a massless photon 
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field. Because of the symmetry between electric and magnetic fields in electrodynamics, 
this system is indistinguishable from one with an electrically charged scalar field, i.e., the 
Abelian Higgs model. 

It is therefore possible that near the critical point, the broken phase of the Georgi- 
Glashow model is dual to the symmetric phase of the Abelian Higgs model. If the duality 
extends through the phase transition, the confining symmetric phase of the Georgi-Glashow 
model is dual to the broken superconducting phase of the Abelian Higgs model, with 
Abrikosov flux tubes playing the role of the confining strings. This is a concrete example 
of the’t Hooft-Mandelstam picture of confinement as a dual phenomenon to superconduc¬ 
tivity I, 0. 

Earlier studies of the Georgi-Glashow model shed some light on this possible dual¬ 
ity 1^^. Best known is the limit A —> oo taken with constant = |m^|/A. It corresponds 
to fixing the norm of the Higgs field <I>, and is traditionally parameterised by couplings 
K = {m? + 8)/A and P = ■ The limits k ^ oo and /3 ^ oo of that theory are, re¬ 

spectively, the compact U(l) gauge theory and the global 0(3) spin model. The former is 


believed to have a weakly first-order phase transition [18|, and the latter a second-order 
one. These two transitions are connected by a phase transition line that separates the Higgs 
and confining phases. There is evidence for a tricritical point at finite (k, /3) at which the 
transition changes from first to second order [31|, but it is not known if the tricritical line 
extends to A = 0. 

Interestingly, the k = oo limit, i.e., the compact U(l) theory, is exactly dual to the 
so-called frozen superconductor . This is an Abelian integer-valued gauge theory, which 
can be obtained as the A —> oo, k —> oo limit of the Abelian Higgs model. In other words, 
the hypothetical duality discussed above is real and exact in this particular limit of the 
theory. The interesting question is whether it exists as an asymptotic duality even away 
from the k = oo limit. This is by no means clear because the compact U(I) theory does 
not have a second-order transition. 

In principle, the methods discussed and used in this paper can be used to test the 
duality hypothesis. If one finds a second-order phase transition, one can measure how the 
masses and the gauge coupling gji change as the transition is approached and determine 
whether the particles decouple. One can then construct observables along the lines 
of Ref. to compare the critical behaviours of the monopoles in the Georgi-Glashow 
model and scalar particles in the Abelian Higgs model. This will, however, be a major 
computation, and is beyond the scope of this paper. 


8. Conclusions 


We have seen how the quantum mechanical mass of a’t Hooft-Polyakov monopole can be 
calculated non-perturbatively using twisted boundary conditions. The method has clear 
advantages over alternative approaches based on creation and annihilation operators and 
fixed boundary conditions. While similar calculations have been carried out before in 
simpler models |W, IT, |T^, this appears to be the first time it has been used for’t Hooft- 


Polyakov monopoles in 3-|-l dimensions. 
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The results demonstrate that one can obtain relatively accurate results for the monopole 
mass. It would be interesting to compare the results with the corresponding classical mass 
to determine the quantum correction. As we have seen, the finite-size effects due to the 
magnetic Coulomb interactions are significant, and therefore one has to compute the clas¬ 
sical mass on the same lattice to have a meaningful comparison. Furthermore, the classical 
mass has to be calculated using the renormalised rather than bare couplings, and this 
introduces a dependence on the renormalisation scheme and scale. 

The simulations in this paper were done at weak coupling, i.e., deep in the broken 
phase. This is a useful limit for testing the method and also for identifying the quantum 
correction. However, the strong-coupling limit, which corresponds to the neighbourhood 
of the transition point, is arguably more interesting. In perturbation theory, the transition 
is of first order, and therefore one cannot reach the critical point, but it is possible that 
this changes if A is high enough. The methods discussed in this paper could then be used 
to study the critical behaviour. 

A particularly interesting possibility is that an asymptotic electric-magnetic duality 
appears near the critical point. The theory would then become equivalent to the Abelian 
Higgs model, with monopoles playing the role of the charged scalars. This would be a 
concrete example of the picture of confinement as a dual phenomenon to superconductivity. 
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